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Abstract. In mesoscopic and nanoscale systems at low temperatures, charge 
carriers are typically not in thermal equilibrium with the surrounding lattice. The 
resulting, non-equilibrium dynamics of electrons has only begun to be explored. 
Experimentally the time-dependence of the electron temperature (deviating from the 
lattice temperature) has been investigated in small metallic islands. Motivated by 
these experiments we investigate theoretically the electronic energy and temperature 
fluctuations in a metallic island in the Coulomb blockade regime, tunnel coupled to 
an electronic reservoir, i.e. a single electron box. We show that electronic quantum 
tunnelling between the island and the reservoir, in the absence of any net charge or 
energy transport, induces fluctuations of the island electron temperature. The full 
distribution of the energy transfer as well as the island temperature is derived within 
the framework of full counting statistics. In particular, the low-frequency temperature 
fluctuations are analysed, fully accounting for charging effects and non-zero reservoir 
temperature. The experimental requirements for measuring the predicted temperature 
fluctuations are discussed. 
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1. Introduction 

A macroscopically large system, thermally coupled to a heat reservoir, is in thermal 
equilibrium with the reservoir when the system temperature is equal to the one of the 
reservoir. Microscopic fluctuations of heat between the system and the reservoir will 
not modify the thermal equilibrium. In smaller, meso or microscopic systems this is not 
necessarily true [1]; individual events of energy transfer between the reservoir and the 
system can induce fluctuations in time of the system temperature [2]. Consequently, 
microscopic heat fluctuations (no external force applied) can drive the small system away 
from thermal equilibrium with the reservoir. The great interest in nanoscale systems 
during the last decades has motivated extensive efforts [31 m 0 El [71E] to reformulate 
or invent novel concepts in order to properly describe the thermodynamics of small 
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systems. Here, as an important example, the concept of temperature fluctuations has 
led to considerable debate [91IIQ]. 

In mesoscopic or nanoscale solid state systems coupled to electronic reservoirs and 
kept at low temperatures, quantum tunnelling of electrons between the system and 
the reservoirs is an important source of individual, random events of energy transfer 
A particularly interesting situation occurs when the electrons in the 
system interact much stronger with each other than with the lattice phonons, often the 
case in both metals |5] and semiconductors [Mj. The electron-electron interactions lead 
to a rapid thermalization of the electron gas, on a time scale Tg-e- As a consequence, 
the electrons are in local thermal equilibrium, characterized by a Fermi distribution 
with temperature Tg. However, the weak electron-phonon interaction leads to that the 
electron gas is thermalized with the lattice phonons on a much slower time scale rg_ph. 
If the energy exchange between the system and the reservoir, due to tunnelling, occurs 
on an intermediate time scale te, i.e. 

Tg.g -C Te -C Tg.ph (1) 

the system electron temperature will develop fluctuations, Tg = Tg(f), with the dynamics 
driven only by quantum tunnelling. 

Electronic temperature fluctuations in mesoscopic systems are typically challenging 
to measure. However, in a number of experiments, focusing on temperature based 
photon detection, dZl [IHl ng 1201 [21] the electron temperature Tflt) of a small metallic 
island was monitored in real time following a pulsed excitation. In the very recent 
experiments in Refs. [221E] , temperature fluctuations of the order of 10 mK away from 
the lattice temperature ~ 100 mK could be detected, approaching the limit of single 
microwave photon detection. Moreover, the observed relaxation time rg_ph of Tflt) 
towards the lattice temperature was of the order of 100 /rs, several orders of magnitude 
longer than the typical rg_g, of the order of 1 ns or below [22] • This large separation of 
time scales puts in prospect experiments with real time monitoring of Tflt), driven by 
electron tunnelling occurring on an intermediate time scale te, fulhlling the inequality 
in Eq. Q. 

In this work we present a theoretical investigation of the tunnelling-induced 
electronic heat transfer and temperature fluctuation statistics in a conceptually simple 
system. Motivated by the experiments HZl ESI nHll2Dl ED we consider a metallic island 
in the Coulomb blockade regime, tunnel coupled to a reservoir, i.e. a single electron 
box [2S|, see Fig. We point out that several theoretical investigations of temperature 
fluctuation distributions have been performed in related systems, e.g. in open, non¬ 
interacting islands [21|, overheated single electron transistors [25], time driven systems 
[26] . superconducting heterostructures IZ71, and an island with injection of electronic 
wave packets [2H]. In our work we focus on the case with no external force applied 
to the system, such as a static or time-varying voltage or thermal bias. In fact, due 
to Coulomb blockade, only zero or one excess electrons are allowed on the island and 
hence, there can be no net charge transfer between the system and the reservoir. These 
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simplifying assumptions allow us to provide a detailed, analytical description of the 
statistical properties of the heat transfer and temperature fluctuations. 

To provide a physically compelling picture, we hrst investigate the statistics of 
energy transfer in thermodynamic equilibrium, when the system electron temperature Te 
is constant and equal to the reservoir temperature Tl. We show how the fluctuations of 
transferred energy depend both on the fluctuations in number of tunnelling events during 
the measurements as well as on the fluctuations in energy transferred in each tunnelling 
event. Based on the understanding of the energy fluctuations in thermal equilibrium, 
we analyse the fluctuations of the system temperature away from equilibrium (Te 7 ^ Tl), 
induced by the fluctuations of the energy transfer. Employing a Boltzmann-Langevin 
approach [291 EQ] for the fluctuation correlations as well as a semiclassical stochastic 
path integral technique [SU |32] for the full distribution of fluctuations, where the 
quantum tunnelling of electrons act as the generator of fluctuations, we analyse the 
statical properties of the temperature fluctuations ATe(f) = Tflt) — {Te(f)). 

For the magnitude of the fluctuation correlations, of direct experimental relevance, 
the result is of the canonical form for fluctuations around thermodynamical equilibrium 

ID 

{ATflt)ATflt + T)) = (2) 

where the correlation time tq = C^/k with the heat capacity of the island and 
K the thermal conductance of the island-reservoir contact, both taken at the reservoir 
temperature Tl. We note that the instantaneous fluctuations are given by [ 1 ] {AT^{t)) = 
k^Tl/Ce, independent of both the island-reservoir contact properties as well as charging 
effects. In contrast, the low-frequency correlator J dt'{ATflt)ATflt + t')) = 2/cbTl/k, in 
accordance with the fluctuation-dissipation theorem for heat transport [231EO], depend 
on both contact properties and charging effects via k. 

Higher-order temperature correlators as well as the full counting statistics, clearly 
demonstrating the non-Gaussian nature of the fluctuations, are investigated with the 
focus on the low-frequency regime. As a general result we hnd that the temperature 
fluctuations increase for increasing charging effects, a consequence of both a wider 
range of energies of tunnelling electrons participating in the heat transfer and larger 
fluctuations in the number of tunnel events during the measurement. 

The rest of the paper is organized as follows. In Sec. 2 the system and model are 
presented. Then, in Sec. 3 the probability distribution of the energy fluctuations are 
derived. In Sec. 4 the energy transfer statistics for equal temperature of the system 
and the reservoir is analysed. Thereafter, in Sec. 5, the distribution of the temperature 
fluctuations is derived and analysed in detail. Finally, in Sec. 6 we conclude and discuss 
prospects for an experimental realization of our proposal. 
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We consider the system shown in Fig. A metallic island, or dot, is conpled via 
a tnnnel barrier, with resistance R and capacitance Cl, to an electrically gronnded 
reservoir. The dot is farther conpled capacitively, with strength Cg, to an electrostatic 
gate, kept at a potential Vg. The dot is in the Conlomb blockade regime, R ^ h/e^, with 
a charging energy E^. = e^/2(CL+Cg). The electronic reservoir is kept at thermodynamic 
eqnilibrinm, at a constant temperatnre Tl, and characterized by a Fermi distribntion 
f{E, Tl) = 1/{l+e'xp[E/k-QTifl) for the electrons, with energy E connted from the Fermi 
energy. 




f(E,p) f(E,TL) 




Figure 1. a) Schematic of the single electron box. A metallic dot is coupled 
capacitively, Cg, to an electrostatic gate, kept at a voltage Vg. The dot is further 
coupled via a tunnel barrier, resistance R and capacitance Cl, to a metallic reservoir. 
The reservoir, in thermodynamic equilibrium, is grounded and kept at the lattice 
temperature R. In the dot, the electrons are in local thermal equilibrium at a 
temperature Tflt) fluctuating in time, b) The Fermi distributions of the electrons 
in the dot (left) and the lead (right), shown schematically. Tunnelling processes in 
and out of the dot with respective rates. Tin and Tout) are shown, c) Representative 
time trace of fluctuation of electron temperature Tflt). Magnification of short time 
interval, with individual tunnelling events visible, shown in rectangle. Typical time 
scales for fluctuation correlations tq and energy tunnelling rg together with average 
temperature (Tg) = R shown. 

On the dot, the electron-electron interaction time re_e is much shorter than the 
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electron-phonon interaction time Te-ph- On time scales longer than Te-e, the electron 
distribntion is thermalized, characterized by a Fermi distribntion f{E, Tg). The electron 
temperatnre Tg = Tflt) generally depends on time and is typically different from the 
constant lattice temperatnre Tl (same as the reservoir temperatnre). The total dot- 
reservoir system is thns in a non-eqnilibrinm state, with the dot electrons in qnasi- 
eqnilibrinm. 

Electrons tnnnel seqnentially in to and ont of the dot with rates and 
T^t respectively (co-tnnnelling is not considered). We consider the sitnation with 
ksTe ^ Ec, and a gate voltage Vg tnned so that only zero or one excess electron 
on the dot is energetically allowed. The time scale rg for energy exchange between the 
dot and the reservoir is given by the typical time for an in-tnnnelling and a snbseqnent 
ont-tnnnelling (or vice versa), i.e. 


_ 1 ^ 1 

in ^ out 


(3) 


For te S> Te-ph the electron-phonon scattering is much faster than the process of energy 
transfer via tunnelling and the dot electron temperature is constant and equal to the 
lattice temperature, Tg = Tl. The dot-reservoir system is thus in thermal equilibrium. 
Below, this regime is only considered in order to provide a qualitative understanding of 
the dot-reservoir energy transfer statistics. The regime of main interest is instead the 
opposite one, described by Eq. ([^, where the dynamics of Tg(t) is driven by electron 
tunnelling. We note that the tunnel rates f|°^ and Fq^^ depend on the difference in 
charging energy for zero and one electrons, further discussed below. 

To provide a physically intuitive picture of the mechanism of the temperature 
fluctuations we consider the following scenario: starting with the dot-reservoir in thermal 
equilibrium, Tg = Tl, an in-and-out tunnelling event take place. The energies of the in 
(out) tunnelling electrons, added to (subtracted from) the total dot electron energy, are 
random quantities with probability distributions determined by the dot and reservoir 
Fermi functions, accounting for the difference in charging energy between zero and one 
excess electron on the dot. As a consequence, after the tunnelling events the total dot 
electron energy has typically changed, with the same probability for an increase as for a 
decrease (as Tg = Tl). Since the dot electron temperature is directly related to the total 
energy, via the heat capacity Cg, and the electron gas is thermalized on the short time 
scale Te-e, the dot temperature Tg changes accordingly. Then, with the next in-and- 
out tunnelling event, Tg changes further and over time the dot electron temperature 
develops fluctuations, as shown in Fig. [T] Importantly, for Tg > Tl the probability 
for a net energy transfer out of the dot during an in-and-out tunnelling event is larger 
than the probability for a net in transfer, and vice versa for Tg < Tl. Hence, there is a 
negative feedback mechanism inherent in the model, characterized by the time scale tq, 
providing relaxation of Tg towards Tl and hence preventing large deviations of Tg away 
from Tl. 
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We consider the dot-reservoir energy transfer statistics a time scale r long enough for a 
large number of tunnel events to occur, r S> te, but short enough for the temperature 
to be taken constant, r <C re- The energy transfer statistics of the dot can be described 
within the framework of a master equation [35l IMl ETJ [3H] , generalized to include energy 
counting helds ^ [32l [SU |38] , 


dt Pi{^,t) 




PoiU) 


(4) 


where Po{^,t) and are the generalized probabilities for zero and one excess 

electron on the dot, respectively (note that Po{0,t) + Pi(0,t) = 1) and the rate matrix 


M(0 = 


rj”’ 


r„.(0 -r 


( 0 ) 


(5) 


The counting held dependent rates to tunnel in to or out of the dot are given by the 
standard single charge tunnelling rates [23] modihed to include counting held factors as 

FtoK) fdEll-f(E,T,)]f(E + A,n)e'i‘^ 


ro „.(0 = 


e 2 p 

1 


dP/(P,Te)[l-/(P + A,TL)]e-*«® 


( 6 ) 


with = rin/out(^ = 0). The energy A = Pc(l) - Pc(0) = Efll - 2nfl) is the 

diherence in charging energy Efln) = Efln — Ugfl for one and zero excess electrons on 
the dot and eug = VgCg the gate induced dot charge, with Cg the gate-dot capacitance. 

We are interested in the full distribution of energy, Pt{£) transferred into (e > 0) 
or out of (e < 0) the dot during the time r. The distribution is conveniently expressed 
in terms of a cumulant generating function (CGF) F{^) as 

^ (^) 

The CGF is given by r times the eigenvalue of the rate matrix in Eq. ([^ which 

goes to zero for —)■ 0, as 


- 2 F( 0 /r = F® + Fi°) 



+ 4Fi,(0Fout(0 


( 8 ) 


We note that since the temperature in the dot and the lead are typically diherent, 
Tl 7 ^ Tefl), we have no simple analytical expression for the tunnel rates in Eq. (§ and 
hence not for F{^). 


3.1. Lowest order cumulants 


The different cumulants of the energy transfer are given by successive differentiation of 
the cumulant generating function with respect to The average energy is 


(£) = 


dP 


d{iO 


= T 


5=0 


p(0) p(l) p(0)p(l) 

p(0) p(0) 

■*' out ■*' in 


(9) 
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Here we have introduced for clarity and later convenience 




p(n) _ 
out 




d{iO 

d^T ont 


«=o 


em 


d{-tO^ 


«=0 




dEEfll-f{E,T,)]f{E + A,n) 

f dEE-f{E,T,)[l-fiE + A,n)]. 


( 10 ) 


We point out that in the steady state the probabilities obtained from Eq. (|^, for = 0, 


are If, = 1 - Pf = r|“V(riS + r|“) 


) and hence (e) = - r^Pf; 


hi) 


Since T; 


( 1 ) 


in/out 


are the energy tunnelling rates, our result is consistent with the dehnition of average 
energy current flowing into the dot as 

(^) 


Je = 


T 


( 11 ) 


The current is plotted as a function of Tg/Ti^ in Fig. for a number of different 
A/lksTf,), clearly showing that the energy current flows from the hotter to the colder 
part of the total system and being zero at Te = Tl, as expected. Moreover, the spectral 
distribution of the energy flow for different temperatures is shown in Fig. 



Figure 2. Left: Probability distribution p{E) (normalized spectral distribution) for 
energy of electrons tunnelling between the dot and the reservoir, for Tg <Ti^,Tg = Ti^ 
and Te > Tl. Solid (dashed) line denotes the distribution for in (out) tunnelling, with 
T > 0 for energy added to the dot. Right: Energy current je (in arbitrary units) as a 
function of Te/Tu, for a = A/(2fcBTL) = 0, 2 and 4. 


The magnitude of the fluctuations of the transferred energy, ((e^)) = {{Aefl) = 
((e — {s)Y) = d?'F/d{i^y‘\^=o is given by the expression 


((Ae)^) 

T 


p{2)p(0) I p(0)p(2) _ pp(l)p(l) 

■*' in out W in out in out 


p(0) p(0) 

out ~P in 


9 I pd) F^^^pd) 

I out-*- in -*- in -*- out 




( 12 ) 


Third and higher order cumulants of the transferred energy are found accordingly. 
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To gain a deeper physical understanding of the dot electron temperature fluctuations, 
it is instructive to hrst consider the case where the dot electrons equilibrate with the 
surrounding lattice on a time scale much faster than the average time te for an in-and- 
out tunnel event. In this case we have a constant Tg = Tl, i.e. the reservoir and the dot 
are in thermal equilibrium during the entire time to of fhe measurement. It is then also 
possible to obtain an analytical expression for the tunnel rates in Eq. and hence 
the CGF in Eq. (|^, allowing for a detailed investigation of the energy transport in the 
system. Making use of the relation/(E-|-a;)[l —/(E)] = lf(E + x) — f(E)]/(l — e^^^^^) 
and the Fourier transform of the Fermi distribution we hnd for the rates 

ri „(0 = rS”>A(0, = r™ a(-0, (13) 


with 


and 


do) _ 1 


A 




h(n = .-^gA/2 27r/cBTLsin(^A/2) 

A sinh(7r/cBTL0 


p(0) _ p(0) A/ZcbTl 
J ■*- out in ’ 


(14) 


(15) 


Introducing the dimensionless ratio a = A/(2 A;bTl) and the dimensionless counting held 
A = we can write Eq. (1^ as (putting r = to) 


F(A) = 


Ct 

e^R sinh a 


cosh a — \ sinh^ a + 


TT sin(Aa) 
asinh(7rA) 


(16) 


Importantly, the analytical form of the CGF as well as the corresponding distribution 
of the energy transfer Ptfls) can be understood in the following way: (i) First, for 
a given measurement there occur a large number N ^ 1 oi in-and-out tunnelling 
events. Since the tunnelling is a quantum process, N will huctuate from measurement 
to measurement. The CGF G{x) fo^' fhc probability distribution P{N) of the number 
of in-and-out events is directly obtained by substituting in Eq. (|^ the off-diagonal 
elements rin/out(0 —t Following the same procedure as above this gives 


G{x) = 


tok-QTE tt 
e^P sinh a . 


cosh a — \/ sinh^ a -b 


(17) 


The square root form of the CGF in Eq. (16) is thus a consequence of the huctuations 
in number of tunnel events A. (ii) Second, in every in-and-out tunnelling event, an 
energy E is hrst transferred into the dot and then an energy E' is transferred out of 
the dot. The probability distribution of the energies E, E' is given from the integrand 
of the tunnel rates in Eq. (|^, as 

/(E + A)-/(E) 


p(E) = 


A 


(18) 


and, taking into account that out-tunnelling for E > 0 removes energy from the dot, 
p(E') = p(—E), as plotted in Fig. Since the transferred energy in each event is 
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independent on the number of events N, in Eq- Q can be written as the 

probability distribution for the energy transferred during N statistically independent 
in-and-out tunnelling events, averaged with respect to the probability distribution for 
N. For the CGF F{X), this implies that we can make use of the cumulant composition 
rule and write 


F(A) = G(x(A)), (19) 

where the CGF x(-^) fhe energy transfer during a single in-and-out event is given 
from the moment generating function 


gix(A) = (^J 


= h(A)h(-A) = 


7rsin(AQ!) 


( 20 ) 


a;sinh(7rA) 

the product of the moment generating functions for energy transfer during in and out 
tunnelling. The function h{\) = h{^ = X/k-^Tifl is given in Eq. (15). From Eqs. (0- 


(20) we see that we obtain the full CGF F{X) in Eq. (16). 


Having analysed the physical origin of F{X) in Eq. (16) we make some comment 


on its analytical structure. First, F{X) is an even function of A, following from that 
sin(AQ;)/sinh(7rA) is even in A. Thus, all odd cumulants of the total energy £ are 
zero. As a consequence, the energy distribution Ftfle) is even in e, around zero (since 
Te = Tl). Second, F{X) is also an even function of a, hence the sign of the charging 
energy difference A is not important for the energy transfer statistics. Third, in contrast 
to the statistics of number of events N with the CGF G{x) Eq. 0 , or the well 
studied charge transfer statistics for a system in the transport state [39l |3H], F{X) is 
not periodic in A. This is a consequence of the fact that the energy transferred in each 
tunnelling event is not quantized. 


fll. Cumulants and probability distribution 


Detailed information about the energy transfer is given by analysing the lowest order 
non-zero cumulants. Since all odd cumulants are zero, the average energy transfer 
{e) = 0, i.e. zero average energy flow, as expected for thermal equilibrium. For the 
second cumulant we have 


((Ae)2) = 


d'^F 


d{iXy 


= tl 


A=0 


2(/cbEl)^ a(7r^ 

' sinh(2Q;) 


( 21 ) 


Making use of the composition relation in Eq. (19) and the fact that odd cumulants of 
x(A) are zero, we can also write 


((Ae)^) 


{kBTLf 


dG 

dflx) 


<Piix) 

x=o d{iXf 


A=0 


( 22 ) 


with 


dG 


d{ix) 


x=o 


/cbElck 

°e^i?sinh(2Q() 


to 

Te 


d?{ix) 2(a^ + TT^) 

<i(iA)2 3 


( 23 ) 
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The energy fluctuations are thus equal to the fluctuations of the energy transferred 
during a single in-out tunnelling event, multiplied by the average number of events 
(N) = to/xE during the measurement. The number of events (N) is exponentially 
suppressed for increasing a = A/2A:b7l, i-e. in the Coulomb blockade regime. In 
contrast, the single event energy fluctuations increase with increasing a, a consequence 
of large A making the energy distribution p{E) in Eq. (18) broader. 

Higher-order cumulants provide information about the deviation from a Gaussian 
distribution. While the odd, third cumulant is zero, the fourth cumulant is nonzero and 
given by 


((£^)) = ((Ae)^) 


(^bTl 


5 cosh^ a 

We can also, similar to the second cumulant, write 


[—a^ -I- vr^ -I- (4a^ -|- bvr^) cosh(2a)] 


(24) 


((e^)) = {kETEY 


fl^G 


d{ixY 


x=o 


d^{ix) 


d{iXy 


A=0 


+ 


dG 


diix) 


dfltx) 


x=o 


dX^ 


A=0 


(25) 


where, in addition to the results in Eq. (23), we have 
d'^G to cosh(2a) d^(w) 


d{ixY 


x=o 


2rE cosm a 


dX* 


4(7r^ — a^) 


A=0 


15 


(26) 


Thus, for {{eY) we see that in addition to the fourth cumulant of the single event energy 
fluctuations multiplied by the average number of events (N), there is an additional 
contribution from the fluctuation in number of events times the square of the single 
event energy fluctuations. Note that total plotted in Fig. has the same sign, 

independent of a. 




Figure 3. Left: Second, ((e^)), and fourth, ((e^)), cumulants as a function of 
a = E.I{flkETE)^ normalized to their respective values at a = 0. Right: Logarithm 
of probability distribution PtoG) for two different a = 0,2, as a function of e‘ = 

Turning to the full distribution Ptflfl)-, Eq. ([^, for small energies e we can (due to 
N ^ 1) perturbatively evaluate the full distribution in the saddle point approximation. 
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To exponential accnracy we have 

\nP,fl8) = -tCe + F{C), 

which gives, to fonrth order in energy 

1 


-ie + dF/d^\^=^* = 0 


In Pto = 


+ 


(27) 


(28) 


2{{s^)) 24 ((£ 2))4 

The logarithm of the distribntion for two different a is plotted in Fig. We note that 
for a ~ 1 and below, the fourth order term oc in Eq. (28) is negligibly small, i.e. the 
distribution is effectively Gaussian. 

To summarize the section, we hnd that for equal temperatures in the dot and the 
lead, Tflt) = Tl, the energy transferred fluctuates from measurement to measurement, 
with only the average over many measurements being zero. The fluctuations are a 
consequence of both the fluctuations in energy transferred during a single in-and-out 
tunnelling event and the fluctuations in the number of events during the measurement. 


5. Temperature fluctuations 

As a result of the total energy fluctuations, it is clear that in the quasi-equilibrium regime 
described by Eq. ([^ the electron temperature Tflt) will develop fluctuations in time. 
The mechanism underlying the fluctuations was discussed above, with a typical time 
trace of the fluctuating temperature shown in Fig. [Tj Our aim is to derive a complete 
description of the statistics of the temperature fluctuations. Formally, considering a 
measurement time to, fhe full distribution of the time-dependent temperature in the 
interval [0,to] can be written |3D] as a functional Fourier transform 

po 


F[Te(t)] = / FC(t)exp 


—i 


dt'T,(t')C{t') + S[C(i)l 


(29) 


where the measure PCit) denotes integration over all possible time-dependent counting 
flelds ({t) and 5'[((’(t)] is the generating functional. Performing n functional derivatives 
of ^[^(t)] with respect to C,{t) we get the n-point (irreducible) temperature correlation 
function in the time domain, as 


nrr., ^ 1 5”>5[C(t)] 

{(r,(i,)...re(i„))) - ^ 


(30) 


C(i)=0 


Importantly, for a typical time te satisfying Eq. ([^, in the range 1 ns to 1 /rs, 
resolving individual tunnelling events in time with existing experimental techniques is 
extremely challenging. Thus, the short time temperature correlators in Eq. (30) are 


difficult to investigate experimentally. Instead, the experimentally accessible quantities 
are the low-frequency correlators, in particular the magnitude of the fluctuations, or the 
temperature noise, given by 


((ATe)2) = - 

to 


io rio 



dtdt'{ATflt)ATflt')), ATflt) = Te - (T,) (31) 


0 JO 
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We note that the temperature noise can be written as {{ATefl) = {{0‘^))/to where have 
introduced the time-integrated temperature 

6= r Tflt)dt. (32) 

Jo 

This motivates us to focus our discussion and present analytical results for the 
fluctuations of 6 . The full distribution of 6 , P{0), can be written, similarly to Eq. 
0. in terms of a CGF S{() as (C is constant in time) 




dCe 


-iCe+siO 


(33) 


From S{() the different cumulants are obtained by successive differentiation, in 
particular the average temperature {TJ) = {0)/to = {l/to)dS{()/d{i()\^=o and the 
fluctuations {{ATefl) = {{0‘^))/to = (l/to)d^*S'(C)/d(*C)^lc=o- 

To provide a physically compelling picture of the temperature fluctuations, below 
we first discuss the average and the time-dependent second-order correlation function 
within a Boltzmann-Langevin approach [29l|30]. Then the generating functional S'[(C(f)], 


and hence the full distribution P[Te(t)], Eq. (29), is derived within the framework of the 
stochastic path integral formalism [mil]. We note that the stochastic path integral 
approach can be understood as a formal extension of the Boltzmann-Langevin approach 
to higher order correlators HU. Moreover, it is equivalent to a fully quantum mechanical 
approach taken in the semiclassical limit. 


5.1. Boltzmann-Langevin approach 

As a starting point, we state the relation between the time-dependent temperature Tflt) 
and the total energy of the dot P(t), given by 


E{t) = Up dee [/(e, Tflt)) - /(e, 0)] = z/p 


TT^klr^it) C'e(Te)Te(t) 


(34) 


6 2 

where z/p the dot density of states at Fermi energy and we introduced C'e(Te) = 
r'p7r^/c|Te/3, the heat capacity of the free electron gas at temperature Tg. We also 
note that the rate of change of the dot energy is, by definition, equal to the energy 
current flowing into the dot 
dE{t) 


dt 


=jm 


(35) 


At steady state, dE(t)/dt = 0, we have j^it) = = 0. The expression for the steady 


state energy current = {s)/to was derived below Eq. (10). The condition of no energy 
current, je = 0, is fulhlled for equal dot and reservoir temperatures, Te = Tp. 

Turning to the fluctuations, following the Langevin scheme, we write the energy 
current and the dot temperature as sums of a constant, time averaged part and a 
fluctuating part, as 

jeit) = je + Ajflt) = Ajflt), Tflt) = Te + ATflt) =Tp + ATflt) (36) 
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The energy current fluctuations can be written as 


+ AT,(i) ^ 


= - «ATe(t) 


(37) 


Te=TL 


where are the ’’bare” fluctuations caused by the stochastic nature of the electron 

tunnelling at the dot-reservoir contact and ATe(t)K is induced by the fluctuating dot 
temperature. Here we introduced k = k(Tl), with /r(Te) = —dji^/dTe, by definition the 
thermal conductance of the dot-reservoir contact at the reservoir temperature Tl. As a 
next step we combine the expressions for dE{t)/dt from Eqs. (34) and (35) and then 
insert the expressions for jflt) and Tflt), from Eq. (36), giving to leading order in the 
fluctuating quantities 

dATflt) 


a 


dt 


= - KAn{t) 


(38) 


where we write for short Ce = Ce{Tfl), the heat capacity of the dot at the reservoir 
temperature Tl. Fourier transforming Eq. (38) and solving for ATflu) we have 

Importantly, at equilibrium the fluctuation-dissipation relation applied to heat transport 
states that, in the semiclassical limit /crTl 3> hu, 

{6jfluj)6je{uj')) = 2kBTlK6{uj + J) (40) 


assuming a white noise spectrum of the energy current fluctuations. We can then write 
the temperature correlator in frequency space 

{AUu)AUu')) = ^^ 6(00 + oo') (41) 

In the time domain this becomes 

{ATflt)ATflt + t')) = (42) 

which is just Eq. (|^. As pointed out above, this result is in the canonical form 
for fluctuations of a physical quantity, in quasi-equilibrium, around thermodynamical 
equilibrium [1]. 


5.2. Stochastic path integral formulation 


For an extension of the Boltzmann-Langevin approach to the full statistical distribution 
of temperature fluctuations, we turn to a stochastic path integral formalism. For 
consistency and completeness of our presentation, we give the derivation of the CGF 
in some detail. As a starting point we consider a time scale r much longer than the 
tunnelling energy relaxation time te but short enough for the temperature to change 
only marginally, r tq. First, starting at time t = 0 with an energy Eo in the dot, the 
probability that the energy in the dot at time r is Ti, i.e. that an energy Eiq = Ei — Eq 
has been transferred, is given by Eq. ([^ 


P{E,o) 


1 


dfoe 


iioiE!i—Eo)+TF(^Q,Eo} 


(43) 
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where is the energy counting field. Extending the reasoning to N time bins, we can 
write the conditional probability that the energy Eiq is transferred through the first 
time bin, E 21 during the second etc., as 

1 f 


P{E^o)P{E2i)....P{Enn-i) = 

I’M-! 


X 


(2vr)^ J ■ 

exp j ^ [-iin{En+l - En) + rF(^„, En)] 




(44) 


, n=0 


The probability that at the end of the N time bins, the total energy E^ — Eq has 
been transferred is obtained by integrating over all intermediate energies En with 
n = 1,N — 1, as 


P{Em) = 


dEi...dEM-iP{EiQ)...P{EMM-i)- 


(45) 


The joint probability to have the electron temperatures Tn = Tflnr) at each respective 
time nr, n = 0,..., iV — 1, can now be written as the conditional probability 


Pn{To^ — J ■■ J dEi..dE]s[_iP{Eifl...P{E]siN_i) 

Af-l 

l[6{Tn-nEn)) 


X 


(46) 


n=0 


where the ^-functions impose the relation between dot energy and temperature, Eq. 
(34), at each time and we introduced T{E) defined from the relation E = C'e(T)T/2. 
Note that the probability distribution P{Eq) for the initial energy Eq will not be of 
importance in the long-time limit we are interested in. Writing the ^-functions as 
Fourier transforms over the variables (^o,..., (n-i we obtain 


-Pv(7o,..., Tn-i) — 


(27r)^ 

with the CGF Sn{Co, given by 

o'S’iv(Cov,C]v-i) _ ^ 


dCo-.dCN-ie 


dEo. .dEjsi- 1 d^Q... d^N- 1 


N-1 

-i T(n.T„ + SN{(o,---,CN-l) 


(47) 


^A7-l 


X exp ^ 


, n=0 


(2vr)^ J ' J 

"l^niEn+l — En) + rF^^ni En) + iTC,nP{En) 


(48) 


Using the fact that the total energy/temperature only changes a small amount in every 
time bin we can now take the continuum limit of Eq. (47), giving the expression in Eq. 
(29) for the full probability distribution P[Te(t)], with the generating functional *S'[(C(t)] 
expressed in terms of a stochastic path integral as 


,5[C(i)] _ 


VE{t)Vat)ew{H[^,E]) 


(49) 


where 


H[^,E] = j\t' (^-^e(^')^^ + E[e(t'),i^(t')] +<(^')r(t')) (50) 
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and 'DE{t),'D^{t) denote integration over all paths in energy and connting field space, 
respectively. 


5.3. Saddle point solution and generating functional expansion 

In the long measnrement time limit considered, with a large nnmber of tnnneling events, 
we can solve the path integral in the saddle point approximation, giving 

S|C(«)| = ffr(«).i5*(()] (61) 

where f,*{t) and E*{t) are the solntions of the saddle point equations 6H[^,E]/6f = 0 
and 6H[^, E]/6E = 0. The saddle point equations are conveniently written in terms of 
and T* = T{E*) as 

dTflt) dF{i*(t),r*{t)) 


2iCe{T* 
-2iCe{r* 


dt df* (t) 

dcit) dF{cit),r*{t)) 


+<(«) 


(62) 

(63) 


dt dT* (t) 

where we note that CflT*) = nF7r^A;|T*(t)/3. 

It is not possible to solve these saddle point equations analytically in the general 
case and hence, we do not have a simple expression for the generating functional ^[(^(t)]. 


However, in order to find the temperature correlation functions in Eq. (30) we can 


expand the generating functional order by order in Cif)- This is done by solving the 
saddle point equations for and T*{t) order by order in ^(t), inserting the solutions 
into 5'[(^(t)] and then performing a functional expansion of 5'[(^(t)] with respect to (fit). 
Here we do not provide any details of this calculation, note however that we expand 

(DO 

= ^o(^) + / dr flit — T)f{T) + ... and similarly for T*it), formally a functional 

— OO 

expansion with time-dependent coefficients fn{t)i%iit). For the first two correlation 
functions we find 


{Tflt)) = 




Kit) 


= Tt, 


C=o 


and 


((Te(t)Te(T))) = 


KS[C] 


Kit)Kit') 




C=0 


2Fl 




(54) 


(55) 


where the correlation time tq = —C'e(TL)/E^T and we have introduced 


Fii=- 


1 d'^F 




FiT — — 


to diif)"^ 

d^F 


€=0,Te=TL 


to d{if)dT 


C=0,Te=TL 


((Ae)^) 

2to^B^L 


(56) 


Here, for F^^ we used the result in Eq. (21) and was evaluated directly 
from Eq. ([^, performing the integrals over the relevant generalized tunnel rates 
rin/out(0 differentiated with respect to f and T. Importantly, using the definition 
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K = —dje/dTe\T^=Ti^ = —F^t we see that the results in Eqs. (55) and (56) are in 
agreement with the results obtained by the Boltzmann-Langevin approach and the 
fluctuation dissipation relation. In addition, the relation F^t = -F^i/iAk^Tl) together 
with Eq. (21) provides us with an explicit expression for the thermal conductance at 
equilibrium, 


K = 


((Aefl) /cbTl a(7r2 + 




3e‘^R sinh(2a) 


(57) 


We note that at a = 0, lifted Coulomb blockade, the thermal conductance is given by 
K = £o7l/( 27?), with Cq = 7r^A;|/(3e^) the Lorenz number. This can be understood a 
Wiedemann-Franz relation |121 03] between the heat conductance k and the effective 
series conductance l/(2i?) for the in-and-out tunnelling events (recall that no net charge 
transport takes place). 

Higher-order correlation functions can be obtained using the same procedure, 
however the expressions become long and difficult to analyse in a transparent way. 
We therefore instead focus on the low-frequency regime where simpler and also 
experimentally more relevant results can be obtained. 


5.3.1. Low-freguency statistics In the long-time, low-frequency limit, we can disregard 
the time-dependence of the energy E(t) and the counting fields ^(t) and C(t) in Eq. 


(49). Consequently, the functionals and path integrals turn into ordinary functions and 


integrals respectively and the probability distribution P[Te{t)] can now be written as a 


to 

distribution of the time-integrated temperature 0 = J dtTflt), as in Eq. (33), with the 
CGF S{9) given from 


= // H(i, E)/to = .(T + F(i, E). 


( 68 ) 


Evaluating the double integral in the saddle point approximation and, similar to the 
discussion above, expanding the resulting CGF order by order in C,, we arrive at two the 
lowest order cumulants. 


{e) = ton, 


= im^) = to 


F2 ■ 


(59) 


as expected from integrating the time-dependent correlation functions in Eqs. (54) and 


(55) over time. Writing out the total result explicitly, using Eqs. (21) and (56), we have 


the low-frequency temperature fluctuations 
((02)) 6e2i?TL sinh(2a) 


ik^Tl 


(60) 


to k-Q a{a‘^ + n) K 

We note that this function is monotonically increasing for increasing a, i.e. stronger 
Goulomb blockade results in larger temperature fluctuations. However, since a large 
a S> 1 also leads to an exponentially large time te between tunnelling events, the 
requirement Te -C Te_ph in Eq. ([^ puts a physical upper boundary on the magnitude 
of the temperature fluctuations. 
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For the third cumulant {{0^)), the skewness of the probability distribution, we can 
proceed in a similar way as for the second one. The resulting cumulant can be written 


m) 


.Ek 

F^rj, 




(61) 


where F^^, F^t are given in Eq. (56) and 


%T - — 


d^F 


F^tt — — 


to d{iifdT 

d^F 


to d{i^)dT^ 


? = 0,Te=TL 


5 = 0,Te=TL 


1 k^Tl 2a [a(a^ + tt^) coth(2a) + tt^ 
to e^R 3 sinh(2a) 

— 1 fcg 2a^ [(a^ + TT^) coth(2a) — a] 


to e‘^R 


3 sinh(2a) 


(62) 


Inserting the expressions for F^^t and F^tt into Eq. (61) we arrive at 


{{6^}) _ _ 108sinh(2a) 

to + 

X [a{a‘^ + TT^) cosh(2a) — (2a^ + vr^) sinh(2a)] . (63) 

The third cumulant is plotted as a function of a in Fig. At a = 0, i.e. lifted 
Coulomb blockade, {{0^))/to = 216e"^F^TL/(/c|7r"^), non-zero and positive. Increasing a 
the cumulant changes sign and increases rapidly in magnitude for increasing a, i.e. for 
strong Coulomb blockade. 




Figure 4. Left: Second, {{9{afl))/{{0{Ofl)), and third, {{0{afl))/{{9{0fl)), cumulant as 
a function of a = A/2 A:b?L) normalized with the value at a = 0. Right: Logarithm of 
probability distribution Pto{9) (solid line) and Gaussian approximation (dashed line) 
for Of = 0 (arb. units). The probability distribution Pto{9) shows clear non-Gaussian 
features with a non-zero, positive skewness. 


To obtain a qualitative picture of the third cumulant we note that Eq. (61) can be 
written 


({«=» 3((Ae)=) 


2\ r 


to 






dn d{{AeY) 


(64) 


dTf, dTg 

with the expression evaluated at = Tp. That is, the fluctuations of induce 
fluctuations of the heat conductance and the magnitude of the fluctuations of the 





























Energy and temperature fluctuations in the single electron box 


18 


transferred energy (on the time scale tq) as dn = {dn/dTflj^Te and = 

{d{{Aefl)/dTe)ATe. These secondary, or cascaded (HI US], fluctnations, are correlated 
with the temperatnre fluctuations and responsible for the third cumulant. 

Instead of considering higher order cumulants individually we plot in Fig. as 
an illustrative example, the full probability distribution Ptfld) in the limit a = 0. The 
probability distribution is evaluated numerically, to exponential accuracy in the saddle 
point approximation. From the plot the non-Gaussian nature of the fluctuations is 
clearly seen: while the distribution eventually goes to zero for 0 = 0 (i.e. minus inhnity 
for the logarithm of PtflO))^ it has a long tail for large 9 which becomes increasingly fat 
for increasing a (not shown), consistent with the large second and large, negative third 
cumulant for large a. 


6. Summary and experimental realization 


To summarize, we have presented a theoretical analysis of the energy and temperature 
fluctuations of a single electron box, a metallic dot in the Coulomb blockade regime 
tunnel coupled to an electronic reservoir. The focus has been on the quasi-equilibrium 
regime, where the dynamics of the dot electron distribution, and hence the temperature, 
is driven by single electron tunnelling. We have investigated the different mechanism 
for the fluctuations of the energy transfer and discussed their effect on the temperature 
fluctuations. 

As a key result, we have presented a prediction for the magnitude of the low- 
frequency temperature fluctuations, in Eq. (60). To provide an estimate of the 


magnitude of the fluctuations we consider the case with lifted Coulomb blockade, a = 0, 
a temperature Tl = 100 mK and take a typical resistance R = 100 kC. Inserted into 
the expression for the fluctuations we have ((0^)) 10® K^s, or equivalently, y/((0^)) ~ 

100 /xK/a/Hz. This is of the same order as the measurement sensitivity reported in the 
most recent experiment [201 EU, suggesting that the predicted temperature fluctuations 
could be observed with existing techniques. 

We also mention that very recent, unpublished, data jlB] show a sensitivity down 
towards ~ 10 /xK/a/Hz. Interestingly, in Ref. 1201 it is claimed that a resolution of this 
magnitude would enable single shot detection of single microwave photons (at 20 GHz). 
Consequently, for the single electron box at large Coulomb blockade, Ajk-Q ~ 1 K, 
where the in-and-out tunnelling process is slow, typically te 3> re_ph, individual electron 
tunnelling events at large energies ~ A could be detected, in single shot, by monitoring 
the dot electron temperature in real time. We conclude by noting that the measurement 
of the dot electron temperature (not discussed in our work) could be carried out in the 
same way as in existing experiments (see e.g. [IS]), by non-invasively tunnel couple the 
dot to a superconducting lead. 
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